REPORT  DOCUMENTATION  PAGE 


Fonn  Approved 
0MB  No.  070U188 


Public  reporting  burden  for  *is  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  corrpleting  and 
reviewing  tiie  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for 
Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503 


1.  AGENCY  USE  ONLY 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 


3.  REPORT  TYPE  AND  DATES  COVERED 


9  Jan  97 


The  Hydrogen  Contribution  To  The  Dark  Matter  In  Draco 


6.  AUTHOR(S) 

Marcus  A.  Boyd 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Arizona 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


96-137 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

DEPARTMENT  OF  THE  AIR  FORCE 

AFIT/CIA 

2950  P  STREET 

WPAFB  OH  45433-7765 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


STRYEMSm  n 


ICS  p'otiic 


12a.  DISTRIBUTION  AVAILABILITY  STATEM 

Unlimited 


13.  ABSTRACT  (Maximum 200  words} 


19910117  192 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


20.  LIMITATION  OF  ABSTRACT 


Standard  Form  298  (Rev.  2-89)  (EG) 
Prescribed  byANSIStd.  239.18 
Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


THE  HYDROGEN  CONTRIBUTION  TO  THE 
DARK  MATTER  IN  DRACO 


by 

Capt  Marcus  A.  Boyd,  USAF 


A  Thesis  Submitted  to  the  Faculty  of  the 
DEPARTMENT  OF  PHYSICS 
In  Partial  Fulfillment  of  the  Requirements 
For  the  Degree  of 
MASTER  OF  SCIENCE 
In  the  Graduate  College 
THE  UNIVERSITY  OF  ARIZONA 


1996 


STATEMENT  BY  AUTHOR 


This  thesis  has  been  submitted  in  partial  ftilfillment  of  requirements  for  an 
advanced  degree  at  the  University  of  Arizona  and  is  deposited  in  the  University  Library  to 
be  made  available  to  borrowers  under  the  rules  of  the  library.  In  addition,  a  copy  is  being 
maintained  at  the  Air  Force  Institute  of  Technology  in  fulfillment  of  my  academic 
agreement  with  the  United  States  Air  Force. 

Brief  quotations  from  this  thesis  are  allowable  without  special  permission, 
provided  that  accurate  acknowledgment  of  the  source  is  made.  Requests  for  permission 
for  extended  quotation  from  or  reproduction  of  the  thesis  manuscript  in  whole  or  in  part 
may  be  granted  by  the  head  of  the  major  department  or  the  Dean  of  the  Graduate  College 
wh^  Jn  Jiis  or  Jier  judgment  4he  proposed  -use  -of  the  Tnateml  4s  m  the  interests  of 
scholarship.  In  all  other  instances,  however,  permission  must  be  obtained  from  the  author 
or  the  Air  Force  Institute  of  Technology. 


SIGNED; 


CUS  A.  BOYD,  CaFt,  USAF 


APPROVAL  BY  THESIS  DIRECTOR 


This  thesis  has  been  approved  on  the  date  shown  below: 


LVIO  MELIA,  Associate  Der 


FULVIC)  MELIA,  Associate  Department  Head 
Department  of  Physics,  University  of  Arizona 


Date 


ACKNOWLEDGMENTS 


First  and  foremost,  I  sincerely  thank  my  thesis  advisor,  Dr.  Fulvio  Melia.  It  has 
been  an  honor  and  a  privilege  to  have  worked  with  and  for  him.  The  knowledge  imparted 
to  me  by  him  has  changed  my  perception  of  the  universe  and  has  sparked  my  interest  in  a 
new  and  excitii^  field  of  science.  His  j)ositive  example  as  a  dihgent  scientist,  caring 
teacher,  and  selfless  mentor  has  left  me  with  a  model  to  follow  the  rest  of  my  life.  Special 
thanks  are  in  order  for  Dr.  Joe  Haller  for  his  expertise  and  prograrns  throughout  the 
course  of  my  research.  Also,  special  thanks  are  in  order  for  Dr.  Robert  Kurucz  for  his 
stellar  model  CD-ROM  and  Dr.  Ed  Olszewski  for  his  expert  advise. 

I  thank  my  wife  Diana,  and  our  daughters  Lindsay  and  Shannon  for  their  patience, 
understanding,  and  unconditional  love  during  this  time  consuming  point  in  my  career. 
Though  these  two  years  in  academia  were  challenging  for  us  all,  we  managed  to  make  it 
with  our  love  for  each  other  still  strong.  I  apologize  for^  the  times  we  did  not  spend 
together.  I  am  in  debt  to  you  aU.  I  also  would  like  to  thank  my  mother  for  teaching  me 
the  importance  of  a  good  education  and  her  constant  pr^ers  during  this  rigorous  period. 

As  for  those  in  the  military,  I  wish  to  thank  the  United  States  Air  Force  Academy 
Physics  D^artment  for  selecting  me  for  the  opportunity  to  obtain  a  Masters  of  Science 
fi'om  the  University  of  Arizona  and  having  ftiith  in  me  to  complete  the  program.  Also, 
many  thanks  to  ROTC  Detachment  20's  staff  and  the  Air  Force  Institute  of  Technology. 
Their  administrative  support,  keeping  me  up  to  date  on  Air  Force  issues  and  sincere 
concerns  for  my  femily  and  my  welfare  is  greatly  appreciated.  I  would  also  like  to  thank 
the  men  and  women  of  the  Air  Force  for  performing  the  daily  duties  I  usually  performed 
when  I  had  a  "real"  Air  Force  job. 

Last,  but  certainly  not  least,  I  thank  God  for  giving  me  the  will  and  strength  to 
accomplish  this  goal  in  my  life. 


ABSTRACT 


High-precision  velocity  measurements  of  faint  individual  stars  in  the  Galactic  dwarf 
spheroidal  galaxy  Draco  indicate  that  the  mass-to-light  ratio  in  this  system  may  be  as 
high  as  ~  100.  Various  possible  forms  of  dark  matter  have  been  proposed  for  this  class  of 
objects,  but  a  more  recent  analysis  based  on  the  Projected  Mass  Estimator  technique  seems 
to  suggest  that  the  non-luminous  mass  is  distributed  and  may  trace  the  stellar  population. 
Here,  we  examine  in  detail  the  possible  contribution  of  H2  to  the  dark  matter  content  of 
Draco  and  conclude  that  most  of  its  mass  may  be  in  the  form  of  10  —  100  Mq  H2  clouds 
if  the  typical  dust  grain  size  is  at  least  ~  90  times  larger  (i.e.,  ^  48  //m)  than  in  our  galaxy. 
This  may  be  consistent  with  the  fact  that  the  dominant  dust  destroying  mechanisms  (e.g., 
supernova  shocks)  that  suppress  the  average  grain  size  in  the  galactic  plane  are  absent  in 
typical  dwarf  Spheroidals.  Thus  the  absence  of  dynamical  processes  (such  as  supernova 
explosions)  that  would  normally  induce  massive  star  formation  (which  is  non-existent  in 
Draco)  may  also  account  for  the  large  grain  size.  The  lack  of  any  observed  extinction  in 
this  system  would  be  consistent  with  the  very  low  areal  covering  fraction  (  0. 1  —  1%) 

of  the  clouds.  We  show  that  the  predicted  H  content  in  Draco  is  then  below  current  21 
cm  detection  limits,  and  that  even  NICMOS  on  HST  would  be  hard  pressed  to  detect  the 
prominent  J  and  K  band  emission  lines  of  excited  H2. 
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1.  Introduction 


Draco  is  one  of  the  low  luminosity  galaxies  in  the  vicinity  of  the  Milky  Way.  The 
eight  Galactic  dwarf  spheroidals  (dSph)  are  among  the  lowest  luminosity  galaxies  known, 
with  absolute  V  magnitudes  of,  e.g.,  -8.9  for  Draco  and  Ursa  Minor,  and  -11.4  for  Leo 
I  (see  Hodge  1971).  Although  grossly  similar  in  stellar  content  to  the  Galactic  globular 
clusters,  their  central  stellar  density  is  typically  about  six  orders  of  magnitude  lower.  In 
addition,  detailed  studies  over  the  past  decade  have  revealed  that  dSph  galaxies  posses  a 
more  diverse  set  of  properties  and  are  comprised  of  more  complex  stellar  populations  than 
the  globular  cluster  analogy  would  predict  (e.g.,  Da  Costa  1992). 

It  is  now  widely  accepted  that  at  least  several  members  of  the  dSph  family,  including 
Draco,  have  very  large  mass-to-light  ratios,  which  are  deduced  from  high-precision  veloc¬ 
ities  of  faint  individual  stars  (e.g.,  Aaronson  &  Olszewski  1987;  Mateo  et  al.  1991).  In 
some  cases,  mass-to-light  ratios  {M/ L)  from  this  core  fitting  range  as  high  as  ~  100,  much 
larger  than  those  found  for  globular  clusters  and  for  most  low-luminosity  ellipticals,  sug¬ 
gesting  that  many  dSph  satellites  of  our  Galaxy  contain  dark  matter.  The  dSphs  are  not 
rotationally  supported  (Paltoglou  &  Freeman  1987;  Mateo  et  al.  1991),  and  their  gas  is 
not  (yet)  detectable  (Knapp  et  al.  1978;  Mould  et  al.  1991),  making  this  tedious  method 
of  measuring  the  velocity  dispersion  the  only  kinematic  technique  for  estimating  the  mass. 
It  is  now  generally  agreed  that  these  kinematic  studies  have  established  the  existence  of 
dark  matter  within  these  objects.  However,  addressing  the  question  of  whether  this  matter 
IS  baryonic  or  nonbaryonic,  and  measuring  its  spatial  distribution,  has  been  substantially 
more  difficult. 

Although  it  is  evident  that  the  “large”  velocity  dispersions  of  ~  10  km  s"^  measured 
for  these  objects  are  real  (as  opposed  to  the  1  km  s“^  value  expected  for  M/L=  2),  it  has 
not  been  clear  where  the  dark  matter  resides.  Both  halo  (Pryor  k  Kormendy  1990;  PK) 
and  central  point-like  (Strobel  &  Lake  1994;  SL)  distributions  have  been  suggested  (see 
§3  below).  Simple  escape- velocity  arguments  show  that  more  than  10^  M©  are  needed 
to  gravitationally  bind  all  of  the  identified  stars  in  the  Draco  system.  While  this  simple 
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physical  argument  is  in  agreement  with  the  analysis  done  from  the  velocity  dispersion  and 
King  models,  it  has  also  been  suggested  that  a  single  central  massive  black  hole  could 
provide  the  mass  needed  to  explain  the  kinematic  results.  Whether  or  not  a  central  cusp 
would  be  expected  to  form  near  these  objects  depends  in  part  on  the  formation  history  of 
the  dSphs.  For  example,  if  the  black  hole  formed  via  accretion  due  to  dynamical  friction, 
the  released  energy  would  puff  up  the  dSph  and  no  cusp  would  now  be  present.  As  such, 
the  absence  of  such  a  distribution  is  not  necessarily  conclusive  evidence  against  the  central 
black  hole  scenario. 

The  possibility  that  neutrinos  might  provide  the  missing  mass  in  Draco  (and  Ursa 
Major)  has  also  been  considered  (Gerhard  &  Spergel  1992),  but  has  been  ruled  out  on 
the  grounds  that  for  a  reasonable  neutrino  mass  (mj,  ~  30  eV)  the  required  core  radius 
(~  10  kpc)  is  substantially  larger  (^  1  kpc)  than  the  size  of  the  dSph  and  the  implied  mass 
(~  4  X  10^^  -Mq)  would  make  the  dSph’s  dynamical  friction  decay  time  in  the  Galactic 
halo  significantly  shorter  than  a  Hubble  time. 

In  this  paper,  we  discuss  in  detail  whether  a  significant  contribution  to  the  dark  matter 
might  be  due  to  the  presence  of  molecular  hydrogen.  An  earlier  search  at  21  cm  (Knapp, 
Kerr  &  Bowers  1978)  has  set  a  tight  upper  limit  of  about  68  M©  on  the  enclosed  mass 
of  atomic  hydrogen  in  Draco,  and  searches  for  the  presence  of  H  II  regions  have  thus  far 
proved  inconclusive.  The  dSph  galaxies  also  appear  to  be  free  of  any  detectable  interstellar 
extinction  due  to  dust.  Counts  of  background  galaxies  indicate  a  very  high  degree  of 
transparency  (Hodge  1971).  Our  purpose  here  is  to  determine,  given  the  physical  conditions 
in  Draco,  how  much  molecular  hydrogen  could  have  thus  far  escaped  detection.  We  begin 
m  §§  2  and  3  by  establishing  the  framework  for  our  calculations,  and  then  develop  the 
physical  state  of  Draco  s  interior  in  §§4  and  5,  wherein  we  consider  the  sources  of  ionizing 
and  molecular  dissociating  radiation,  and  the  chemical  network  responsible  for  maintaining 
the  equilibrium  populations  of  the  various  species.  We  describe  our  calculational  procedure 
and  present  the  results  in  §6,  and  then  discuss  the  observational  consequences  in  §  7,  with 
concluding  remarks  in  §8. 
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2.  General  Considerations 


Globular  clusters  (to  which  the  dSphs  have  sometimes  been  compared)  appear  to  be 
devoid  of  interstellar  matter,  primarily  due  to  the  stripping  of  gas  during  their  passage 
through  the  Galactic  plane  (occurring  roughly  every  10®  years  or  so),  and  the  continuous 
outflow  of  the  hot  interstellar  matter  as  it  is  produced  (Scott  &  Rose  1975).  In  the  case  of 
dSph  galaxies,  however,  the  first  mechanism  does  not  apply  and  the  second  relies  on  the 
gas  temperature  which  may  be  very  low.  It  is  not  unreasonable  to  expect  that  dSphs  may 
have  retained  substantial  quantities  of  primordial  gas. 

Recently,  Elmgreen  (1993)  has  demonstrated  that  the  transition  from  H  to  H2  in 
the  interstellar  medium  generally  depends  sensitively  on  the  pressure  and  radiation  field, 
such  that  even  a  slight  increase  in  pressure  can  convert  a  whole  cloud  population  into 
H2  molecules,  whereas  a  slight  increase  in  radiation  can  convert  the  clouds  back  to  H. 
Almost  all  of  the  H2  forms  on  the  surface  of  dust  grains,  whereas  its  destruction  is  the 
result  of  radiative  and  collisional  dissociation  (Shull  &  Beckwith  1982).  This  is  supported 
observationally  by  absorption  line  studies  of  H2  regions  along  the  line  of  sight  to  many 
OB  stars  (Rogerson  et  al.  1973),  which  have  demonstrated  that  the  fractional  abundance 
of  H2  is  determined  primarily  by  the  balance  of  dust  versus  UV  radiation  in  the  ambient 
medium.  Molecular  clouds  form  when  the  cloud  column  density  is  sufficiently  high  that 
the  absorption  lines  become  self-shielding. 

The  key  ingredients  that  feature  in  our  considerations  are  therefore  (1)  how  the  dark 
matter  is  distributed  within  Draco,  i.e.,  what  its  radial  profile  and  volume  filling  factor  are; 
(2)  the  spectrum  of  the  UV  radiation  field,  which  here  will  be  a  composite  of  components 
from  within  Draco,  from  the  extragalactic  background,  and  more  importantly,  from  the 
Milky  Way  Galaxy  itself;  (3)  the  dust  content;  and  (4)  the  detectable  spectral  signature  of 
H  and  H2.  We  will  establish  what  the  most  likely  radial  distribution  of  the  missing  mass 
is  in  the  following  section,  using  the  projected  mass  estimator  (Haller  k  Melia  1996).  The 
volume  Ailing  factor  must  be  treated  as  a  free  parameter.  The  issue  here  is  whether  or  not 
the  gas  condenses  into  clouds,  and  what  their  size  is.  In  this  regard,  the  Jeans  criterion  for 
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collapse  will  give  some  indication  of  the  permitted  range.  Another  indicator  is  provided 
by  the  required  degree  of  shielding  from  the  external  UV  radiation  for  the  cloud  interior. 
Since  the  ratio  of  cloud  surface  area  to  its  volume  scales  inversely  as  the  radius,  the  ratio 
of  H  to  H2  abundance  will  depend  to  some  extent  on  the  cloud  size. 

Establishing  the  dust  content  in  Draco  is  much  more  difficult,  in  part  because  the 
conditions  in  a  dSph  are  very  different  from  those  in  the  galactic  interstellar  medium, 
where  we  are  guided  by  observations  in  inferring  an  average  dust  to  gas  particle  number 
ratio  of  about  (e.g.,  Black  &  Dalgarno  1977)  and  a  typical  grain  size  of  ~  0.5  fxm. 

The  lack  of  any  detectable  dust  (see  §  1  above)  is  a  poor  guide  if  the  areal  filling  factor  is 
small.  Within  our  galaxy,  it  appears  that  as  much  ~  0.3  -  0.5  of  the  C,  N,  and  0  atoms,  in 
addition  to  most  of  the  heavier  elements,  are  concentrated  within  the  grains  (Spitzer  1978). 
However,  given  the  rapid  rate  of  growth  of  the  grain  size  (~  2  x  10~^^  cm  yr“^),  it  is  not  yet 
fully  understood  why  grains  typically  do  not  grow  very  much  beyond  ~  0.5  fim  in  radius. 
It  is  likely  that  the  mechanisms  that  destroy  the  grains  are  sensitive  to  their  size.  Several 
important  examples  of  these  are  sputtering  and  grain-grain  collisions  in  interstellar  shocks, 
sputtering  and  sublimation  in  H  II  regions,  photodesorption  by  UV,  sputtering  by  cosmic 
rays,  and  sublimation  during  a  supernova  radiation  pulse,  which  is  particularly  important 
for  the  more  volatile  mantle  materials.  Photodesorption  and  sputtering  in  shock  waves  are 
the  principal  destruction  mechanisms  for  less  volatile  mantle  ices.  Sputtering  and  grain- 
grain  collisions  account  for  most  refractory  grain  destruction  (Draine  &  Salpeter  1979). 
Interstellar  shocks  are  presumably  much  less  prevalent  in  Draco  than  in  typical  regions 
of  the  galaxy,  mainly  due  to  the  absence  of  supernovae  and  the  (expected)  extremely  low 
rate  of  cloud-cloud  collisions.  In  addition,  the  cosmic  ray  intensity  at  the  location  of  Draco 
~  67  kpc  from  the  galactic  center  is  expected  to  be  significantly  lower  than  it  is  within 
the  galaxy.  This  lack  of  any  important  destruction  processes  may  lead  to  a  growth  of  dust 
grains  beyond  the  level  seen  in  typical  locations  within  our  galaxy.  We  shall  therefore  treat 
the  average  grain  size  Vg  in  Draco  as  a  parameter  to  be  determined,  and  consider  a  range 
of  reasonable  values  up  to  ~  48  fj,m. 
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Closely  coupled  to  this  is  the  metal  abundance  in  Draco.  Using  the  spectral  flux 
distributions  of  red  giants  in  this  system  over  the  range  AA  3240  -  7620,  Zinn  (1978) 
concluded  that  the  mean  metal  abundance  [Fe/H]  =  -1.86  ±  0.09  is  actually  larger  than 
that  m  some  of  the  globular  clusters,  such  as  M92.  We  shall  adopt  this  value  of  the 
metallicity  in  our  calculations  and  we  will  assume  that,  just  like  the  situation  in  our 

galaxy,  the  dust  grains  in  Draco  contain  about  half  of  the  C,  N,  and  0  atoms  and  most  of 
the  heavier  elements. 

Ultimately,  the  amount  of  H  and  Hj  gas  must  be  reconciled  with  the  available  obser¬ 
vations.  In  the  interstellar  medium,  UV  absorption  line  studies  of  H2  regions  are  made 
possible  by  the  nearby  10  -  30  pc)  presence  of  OB  stars,  which  act  as  background 
sources  (Jura  1975).  Dwarf  Spheroidal  galaxies  are  generally  devoid  of  such  UV  bright 
stars,  making  this  method  of  detection  unlikely.  Under  normal  circumstances  in  inter¬ 
stellar  clouds,  the  temperatures  are  not  high  enough  to  excite  molecules  to  vibrational 
or  rotational  states  above  ground.  However,  a  substantial  fraction  of  the  molecules  may 
be  excited  into  vibrational  or  high  rotational  states  by  UV  radiation  (Gould  &  Harwit 
1963),  which  then  de-excite  by  emitting  infrared  radiation.  In  §  7  below,  we  will  calculate 
the  expected  H2  spectrum  of  Draco  and  discuss  the  instrument  sensitivity  required  for  a 

detection.  The  H  content  is  already  limited  to  less  than  about  68  M©  by  the  available  21 
cm  data  (see  above). 


3.  The  Dark  Matter  Component  in  Draco 

a)  Inferring  Mass  Distributions  Using  The  Projected  Mass  Estimator 

One  can  exploit  a  “tracer  population”  with  observed  radial  velocities,  and  pro¬ 
jected  positions,  (xi,i/i),  associated  with  a  gravitating  system  to  infer  its  total  mass  using 
the  projected  mass  estimator  (PME)  (Bahcall  &  Tremaine  1981  BaT;  Heisler,  Tremaine 
&  Bahcall  1985  HTB).  The  average  of  the  quantity  r^/G  over  the  whole  tracer 
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population  is  related  to  the  total  mass  M  by 

^  JL  ^  _  /  32 

^  ^  hi  ^  ~ 

where  the  z-axis  is  oriented  along  the  line-of-sigh 
from  the  center  of  the  distribution  as  projected  onto  the  plane  of  the  sky.  The  normalization 
factor,  C,  depends  on  the  shape  of  the  gravitational  potential  and  the  distribution  of  the 
orbits  of  the  tracer  population.  For  a  Keplerian  potential  (i.e.,  a  central  concentrated 
mass),  It  is  equal  to  1,  2,  or  3,  when  the  orbital  distributions  of  the  tracer  population 
are  respectively  radial,  isotropic,  or  circular.  On  the  other  hand,  when  the  system  is 
self-gravitating  (no  point  mass)  the  value  of  C  is  0.5,  1,  or  1.5  for  the  same  cases  of 
the  distribution  of  orbits.  Given  N  sample  objects  drawn  from  the  tracer  population, 
the  fractional  standard  deviation  of  the  inferred  mass  M  is  proportional  to  N-^^.  This 


)  r^)  =  CM  ,  (1) 

it  and  r_Lt  is  the  radius  of  tracer  object  i 


procedure  has  certain  advantages  over  the  virial  theorem  (BaT),  yet  one  limitation  is  that 

proper  motion  data  are  generally  unavailable  and  assumptions  are  necessary  to  constrain 
the  value  of  C . 

Haller  and  Melia  (1996)  showed  how  one  could  use  the  projected  mass  estimator  to  infer 

not  only  the  total  mass,  but  also  the  mass  distrihuUon  of  a  gravitating  system.  Consider 

a  spherically  symmetric  distribution  of  tracer  objects  in  steady  state  with  no  rotation 

and  a  known  distribution  function,  /,(r,v),  and  allow  for  the  presence  of  a  spherically 

symmetric  distribution  of  matter  due  to  other  components,  M2(r).  The  velocity  dispersion 

ellipsoid  for  the  tracer  objects  will  be  aligned  with  the  radial  coordinate  basis  vectors. 

Since  the  system  is  spherically  symmetric,  the  transverse  velocity  dispersions  are  equal: 

<re\r)  =  a^^{r)  =  ^^^(r).  Let  us  define  the  quantity  A(r)  =  at^{r)/ar\r\  which 

describes  the  shape  of  the  velocity  ellipsoid  as  a  function  of  radius.  For  a  given  quantity  X, 

the  expectation  value  <  AT  >  is  weighted  by  the  volume  density  of  the  tracer  population: 

<  X  >  =  JyXpidr/  fypidr.  Thus,  the  quantity  {v^^rj_)  for  an  arbitrarily  specified 
volume  V  is  given  by 

{v,^rj_)v  =  ^  r  sin^pi(r)(T^^(r)[cos^  6  +  A(r)sin2  9]  dr 

/y /i(r,v)drdv  ~  /vPi(r)dr  ' 
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Using  the  above  relation  and  the  spherical  Jeans  equation,  it  can  be  shown  that  tor  a 

sample  of  tracer  objects  contained  within  a  circular  aperture  centered  on  the  distribution, 
a  more  general  form  of  Eq.  (1)  is 


The  averaged  quantity  on  the  left  hand  side  of  Eq.  (3)  can  be  readily  obtained  from  obser¬ 
vations.  If  one  has  an  underlying  model  of  the  distribution  function  of  the  tracer  objects 
and  the  enclosed  mass  function  of  other  components  MaCr),  then  all  of  the  quantities  on 
the  right  hand  side  can  be  computed.  Each  of  the  terms  on  the  right  is  a  ratio  of  integrals 
and  are  functions  of  the  projected  radius  R  of  the  aperture  that  contains  the  volume  V  of 
tracer  objects.  We  call  such  a  sampling  the  “aperture  PME”.  In  the  limit  as  R oo  one 


recovers  the  relation  given  by  Eq.  (1). 

We  call  the  factor  C  the  orbit-type  ’  factor  because  it  is  an  average  of  the  velocity 
ellipsoid  shape  parameter  A(r),  weighted  by  the  radial  kinetic  energy  density  of  the  tracer 
population.  ksR-^oc,  it  generalizes  the  factor  C  appearing  in  Eq.  (1).  The  first  term 
inside  the  brackets  of  Eq.  (3)  is  a  density  weighted  average  of  the  enclosed  mass  function 
for  the  tracer  population  taken  over  the  sampled  volume  (in  this  case  a  cylindrical  volume). 
Similarly,  the  term  <  M2(r)  >  is  a  density  weighted  average  of  other  mass  components  but 
the  weighting  comes  from  the  tracer  population,  not  the  density  of  these  other  components. 


Finally,  the  last  term  we  call  the  “excess  measured”  mass  Mb.  It  arises  because  one  is 
sampling  only  a  finite  volume  and  this  boundary  term  is  nonvanishing  at  finite  radii.  It 
may  even  be  quite  important  at  smaller  radii.  When  the  sampling  volume  is  a  projected 
cylinder,  the  variation  of  these  terms  with  R  is  quite  different  from  their  variation  with 
radial  coordinate  r  if  one  were  to  use  a  spherical  sampling  volume. 

The  standard  deviation  of  the  projected  mass  estimator  for  the  given  sample  volume 

IS  given  by  iV  /  {Z2/nG)^j\v:r‘^r\)v  -  where  N  is  the  number  of  sampled 

objects.  (It  is  assumed  that  the  distribution  of  sampled  objects  reflects  the  intrinsic  dis¬ 
tribution.) 
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b)  Application  Of  The  Aperture  PME  To  Draco 

The  dark  matter  distribution  for  Draco  was  modelled  by  PK  using  two-component 
King  models  (1966)  to  estimate  the  mass  of  a  possible  extended  halo.  An  alternative 
model  for  this  galaxy  was  proposed  by  SL,  who  assumed  that  the  gravitational  potential 
IS  dominated  by  a  central  black  hole  with  a  mass  of  order  10^  Mq.  These  two  models 
represent,  in  a  way,  extreme  cases  for  the  mass  distribution  for  this  galaxy.  The  best 
fitting  suite  of  models  presented  by  PK  show  that  it  is  possible  for  the  dark  matter  to 
be  more  extended  than  the  visible  population  of  stars.  For  simplicity  of  comparison,  we 
restrict  our  attention  to  the  PK  model  which  assumes  the  dark  matter  is  distributed  in 
the  same  way  as  the  stars  and  the  total  ratio  of  dark  to  light  matter  is  100.  This  same 
M/L  ratio  is  also  assumed  for  the  model  proposed  by  SL. 

For  a  King  model,  the  gravitational  potential  and  energy  per  unit  mass  E  are 
customarily  measured  relative  to  the  tidal  radius  r*  of  the  system: 

=  -D  +  t/o  and  e  =  -E  +  Uq  for  Uq  =  (4) 

The  distribution  function  of  the  tracer  stars  is  taken  to  be  a  truncated  Maxwellian  given 

by 

/i(e)  =  (27r(7|-)-t  -  1)  e  >  0,  (5) 

while  /i(e)  -  0  for  e  <  0.  For  units  such  that  the  scale  radius  tq,  G,  and  central  density, 
Po,  are  unity,  then  the  square  of  the  scale  velocity  <7^  has  a  dimensionless  value  of  47r/9. 

The  local  density  is  a  function  of  the  local  potential  via  an  integration  of  Eq.  (5) 
over  all  velocities  up  to  the  escape  velocity,  The  one-dimensional  velocity  dispersion 

a^(r)  is  also  computable  (see  BT,  p.235).  Choosing  a  ratio  and  setting  its  initial 

gradient  dxf/dr  =  0,  Poisson’s  equation  is  integrated  to  obtain  the  functions  xf{r),  pi(r), 
and  Mi(r).  The  King  model  parameters  obtained  by  PK  to  fit  the  Draco  system  are  ^o/^- 
—  3.15,  scale  radius  tq  =  156  pc,  and  scale  velocity  =  11.4  km  sec~^. 

Strobel  and  Lake  (1994)  proposed  a  distribution  function  of  tracer  stars  that  resulted 
from  the  formation  of  a  massive  black  hole  of  order  10"  Mq  from  a  primordial  isothermal 
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distribution  of  stars  with  uniform  density  pQ  and  one- dimensional  velocity  dispersion  cxy. 
As  the  black  hole  formed,  the  stellar  distribution  function  was  an  adiabatic  invariant 
(Young  1980)  and  those  stars  that  remained  bound  to  it  had  initial  energies  Ei  << 
which  results  in  the  distribution  function 


m 


Po 

(27r(7„2)3/2- 


(6) 


SL  set  (7i;  =  3  km  sec  ^  based  on  the  typical  velocity  dispersion  of  Galactic  globular 
clusters.  The  density  of  stars  is  then  found  by  integrating  the  distribution  function  in  Eq. 
(6)  over  all  velocities  up  to  the  escape  velocity  at  each  radius,  while  the  velocity  dispersion 
(assumed  to  be  isotropic)  (J|;^(r)  is  also  an  analytic  function.  By  fitting  the  projected 
surface  density  and  velocity  dispersion  profiles,  they  obtained  a  tidal  radius  for  Draco  of 
Tt  =  720  pc  and  a  central  black  hole  mass  M.  (1.5  -  2.0)  X  10^  Mq. 

To  test  which  of  these  models  is  better  able  to  account  for  the  stellar  velocity  data, 
we  computed  the  aperture  P ME  as  a  function  of  R  using  the  best  fitting  parameter  values 
given  by  these  workers.  The  resulting  aperture  PME  curves  were  then  compared  with  the 
radial  velocity  data  set  of  Olszewski  (1994)  (Figure  1).  Although  the  statistical  errors  are 
large,  the  case  for  an  extended  dark  matter  distribution  in  Draco  seems  more  favorable. 
In  particular,  if  a  central  point  mass  were  dominating  the  gravitational  potential  for  this 
system  then  one  would  not  expect  the  PME  to  decrease  so  strongly  hs  R  0.  Generic 
models  that  incorporate  a  massive  black  hole  as  the  dominant  mass  component  all  have  a 
high  value  for  the  aperture  PME  when  the  radius  of  the  sample  aperture  is  less  than  the 
size  of  the  core  radius  of  the  tracer  population  (Haller  &  Melia  1996).  On  the  basis  of  these 
results,  we  conclude  that  the  dark  matter  in  Draco  is  most  likely  distributed,  and  may  in 
fact  trace  the  luminous  matter.  In  Figures  2  and  3,  we  show,  respectively,  the  dark  mass 
density  distribution  for  the  best  fit  PK  model  (see  above)  and  the  corresponding  scales  of 
velocity. 


4.  The  Radiation  Field 

The  radiation  field  within  Draco  is  produced  primarily  by  (i)  its  stellar  content,  (ii) 
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the  extragalactic  background  and  (iii)  the  stars  within  our  own  galaxy.  Although  these 
sources  emit  radiation  across  a  broad  spectral  range,  the  primary  band  of  radiation  that 
severely  affects  the  detectability,  the  structure  and  the  composition  of  a  typical  gas  cloud 
is  the  far-ultraviolet  band  (FUV).  This  band  has  a  wavelength  range  of  900-3000  1  and  is 
responsible  for  the  excitation,  the  photodissociation  and  ionization  of  important  molecular 
species.  Also,  it  is  responsible  for  the  excitation  and  ionization  of  some  atomic  species. 
Specifically,  when  hydrogen  in  atomic  (H  I)  and  molecular  (H2)  form  are  the  dominant 
species  (which  is  usually  the  case  for  galactic  gas  clouds  with  a  cosmic  abundance),  it  is  the 
Lyman  and  Werner  bands  (912  -  1108  A)  that  are  the  most  important  (Shull  &  Beckwith 
1982;  Black  &  VanDishoeck  1987;  Martin,  Hurwitz  &  Bowyer  1990;  Sternberg  &  Dalgarno 
1989)  since  H2  is  primarily  excited  and  dissociated  to  form  H  I  by  FUV  photons  in  these 
bands.  Excitation  occurs  initially  by  a  Lyman  or  Werner  band  photon  exciting  H2  from  its 
ground  electron  state,  ,  to  a  higher  electron  state  of  (for  Lyman)  or  (for 

Werner).  This  stage  of  the  excitation  process  is  followed  by  a  fluorescence  90%  of  the  time 
to  a  bound  vibrational-rotational  excited  ground  electron  state  of  H2.  For  the  remaining 
10%  of  the  time,  dissociation  occurs  by  H2  fluorescence  to  a  vibrational  continuum  ground 
state  to  form  H  I  (Black  &  VanDishoeck  1987;  Martin,  Hurwitz  &  Bowyer  1990). 

The  FUV  radiation  field  due  to  the  stellar  content  of  Draco  itself  may  be  safely  as¬ 
sumed  to  be  negligible  compared  to  the  other  sources  of  FUV  photons  impingent  on  Draco 
from  outside  the  dwarf  Spheroidal  Galaxy.  The  primary  argument  for  this  is  that  the 
stellar  content  of  Draco  is  heavily  biased  away  from  stars  producing  FUV  radiation.  FUV 
producing  stars  are  primarily  those  with  a  MK  spectral  type  of  O,  B  or  A.  Observations 
of  dwarf  Spheroidals  (dSph)  indicate  that  these  types  of  stars  are  extremely  rare  in  these 
systems  and  might  only  exist  on  the  blue  portion  of  the  horizontal  branch  of  a  dSph’s  color 
magnitude  diagram  (Olszewski  1994).  This  is  clearly  evident  on  the  horizontal  branch  of 
the  color  magnitude  diagram  constructed  for  Draco  by  Stetson  (1980).  From  this  portion 
of  the  diagram,  only  2  out  of  365  stars  used  to  produce  the  diagram  may  be  considered 
early  to  late  A  type  and  none  are  O  or  B  type  stars.  The  only  candidate  for  an  0  or  B 
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type  star  on  the  entire  diagram  is  a  peculiar  star  off  the  principal  sequence,  but  it  is  at 
most  an  extreme  late  B  type. 

The  extragalactic  background  is  primarily  produced  by  distant  galaxies,  quasars  and 
the  intergalactic  medium  (Paresce  &  Jakobsen  1980).  The  extragalactic  FUV  background 
has  been  investigated  in  detail  by  Martin,  Hurwitz  &  Bowyer  (1991),  who  report  an  upper 
limit  to  its  intensity  of  280±35  photons  cm"^  sr“^  within  a  band  of  1400-1900  A. 

Assuming  the  background  intensity  to  be  isotropic  within  Draco  and  conservatively  taking 
the  upper  limit,  we  estimate  a  total  flux  of  ^  3.6  x  10“^  photons  cm“2  Hz"^  within 
this  band.  As  we  shall  see  below,  this  component  of  the  FUV  radiation  within  Draco  is 
also  negligible  compared  to  the  contribution  from  the  Milky  Way,  and  it  may  therefore 
itself  be  safely  ignored. 

To  calculate  the  flux  of  FUV  photons  from  our  galaxy,  we  have  developed  a  simple,  yet 
fairly  accurate,  model  of  the  galactic  sources  to  produce  the  spectrum  observed  at  Draco 
(Figure  4).  This  model  consists  of  the  various  spectral  types  and  luminosity  classes  of  stars 
arranged  uniformly  on  a  thin  plane  which  we  call  the  stellar  disk.  The  stellar  spectral  types 
range  from  05  to  M5  and  the  luminosity  classes  span  the  range  from  the  main  sequence  to 
the  super  giant  branch.  The  radius  of  the  stellar  disk  is  taken  to  be  15  kpc  (Allen  1973). 
For  the  purpose  of  incorporating  the  effects  of  line-of-sight  attenuation,  we  surround  the 
stellar  disk  with  a  thicker  disk  of  dust,  which  is  also  distributed  uniformly.  Its  thickness  is 
taken  to  be  300  pc,  in  agreement  with  the  effective  thickness  of  our  galaxy  for  interstellar 
absorption.  Relative  to  this  galactic  plane,  Draco  lies  at  a  heliocentric  distance  of  67  kpc 
(Hodge  1971),  a  galactic  latitude  of  35  degrees,  and  a  galactic  longitude  of  86  degrees. 
Assuming  a  distance  of  8.5  kpc  to  the  Galactic  center,  this  yields  about  the  same  distance 
(fs  67  kpc)  from  Draco  to  the  center  of  our  galaxy. 

The  integrated  galactic  spectrum  is  assembled  from  the  spectra  of  individual  types  of 
stars.  This  is  done  by  binning  the  stars  within  our  galaxy  into  a  MK  spectral  type  of  05, 
BO,  B5,  AO,  A5,  FO,  F5,  GO,  G5,  KO,  K5,  MO  or  M5  within  a  luminosity  class  of  super 
giant,  giant  or  main  sequence.  For  simplicity,  white  dwarfs  and  sub  dwarfs  are  classified 
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as  main  sequence  stars,  and  sub-giants  are  assumed  to  be  giants  and  bright  giants  are 
assumed  to  be  super  giants.  Based  on  the  available  data  (e.g.,  Mihalas  &  Binney  1981; 
Allen  1973),  an  average  effective  temperature  (Teff),  surface  gravity  (g)  and  radius  (Rs) 
may  be  assigned  to  each  spectral  type.  With  these,  and  the  assumption  of  solar  metallicity, 
we  can  then  choose  among  the  Kurucz  (1993)  stellar  atmosphere  model  spectra  for  each 
stellar  type,  taking  the  closest  match  between  Teff  and  g.  Tables  lA,  IB,  and  1C  show 
the  different  parameters  used  for  each  model  and  their  originally  assigned  values. 

To  determine  the  number  of  stars  within  each  type,  we  assume  a  total  galactic  stel¬ 
lar  content  of  ~  1  x  10^^  stars  and  then  use  the  luminosity  functions  in  Bahcall’s  (1986) 
galaxy  model  to  determine  the  number  of  stars  at  each  absolute  visual  magnitude.  These 
luminosity  functions  consist  of  a  disk  and  spheroidal  distributions  with  an  absolute  visual 
magnitude  ranging  from  -7  to  16  in  each  case.  The  use  of  these  luminosity  functions  has 
been  validated  recently  by  the  demonstration  that  they  agree  very  well  with  observational 
data  taken  with  the  Hubble  telescope  (Bahcall  &  Maoz  1993).  The  total  luminosity  func¬ 
tion  for  the  stellar  disk  is  taken  to  be  the  sum  of  these  two  distributions  at  each  of  the 
corresponding  magnitudes.  Merging  this  total  luminosity  function  with  a  spectral  type 
distribution  (e.g.,  Allen  1973,  corrected  by  Bahcall  &  Soneira  1981)  covering  an  absolute 
visual  magnitude  range  between  -7  and  16,  it  is  straightforward  to  derive  a  spectral  lumi¬ 
nosity  function  from  which  we  infer  the  percentages  of  each  spectral  type.  The  different 
spectral  types  are  then  categorized  according  to  their  visual  magnitude  using  data  pro¬ 
vided  for  each  luminosity  class  in  Allen  (1973),  and  Mihalas  &  Binney  (1981).  From  this 
categorization,  and  assuming  that  the  percentages  increase  logarithmically  from  05  to  B5, 
and  adjusting  the  percentages  of  AO  and  A5  stars  to  match  the  HR  diagram  of  the  solar 
neighborhood,  we  infer  the  percentage  of  each  spectral  type  for  a  given  luminosity  class  as 
shown  in  Figure  5  (see  also  Tables  2A,  2B,  and  2C). 

The  attenuation  due  to  dust  and  gas  is  calculated  by  using  the  average  interstellar 
attenuation  curve  of  Savage  &  Mathis  (1979),  which  is  shown  in  Figure  6.  Here,  Rx  is  the 
ratio  of  the  magnitude  of  attenuation  Ax  at  a  certain  wavelength  to  E(B  —  V).  We  note. 
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however,  that  the  dust  within  the  galaxy  is  patchy,  and  one  can  at  best  therefore  only 
determine  the  average  value  of  E{B  -  V\  and  hence  Ax.  We  do  this  by  determining  the 
average  distance  traversed  by  a  path  through  the  disk  of  dust  from  sources  in  the  stellar 
disk  toward  Draco.  Then,  by  assuming  a  visual  attenuation  rate  of  1.9  mag/kpc  (Allen 
1973)  and  of  3.2  (Mihalas  &  Binney  1981),  we  infer  an  average  value  of  E{B  -  V) 
equal  to  0.156  by  multiplying  the  average  distance  with  the  visual  attenuation  rate  and 
dividing  by  R^.  This  average  E{B  —  V)  yields  the  magnitude  of  the  attenuation  at  each 
given  wavelength.  We  also  assume  that  all  the  radiation  with  a  wavelength  less  than  912  X 
(i.e.,  the  Lyman  limit)  does  not  escape  from  the  galaxy. 

The  total  FUV  spectrum  at  Draco  is  the  sum  of  the  spectrum  from  each  star  at  a 
distance  of  67  kpc  from  the  stellar  disk,  attenuated  by  the  factor  derived  above.  The 
result  is  shown  in  Figures  7  and  8  in  units  of  photons  cm-2  s"!  Hz"!  and  photons  cm-^ 

A  ,  respectively.  It  is  evident  that  the  FUV  flux  produced  at  Draco  by  the  Milky 
Way  easily  dominates  over  the  contributions  from  the  extragalactic  background  and  stars 
within  Draco  itself.  In  Figure  9,  we  compare  our  spectrum  with  those  calculated  and 
observed  for  the  interstellar  medium,  which  are  expected  to  show  a  similarity  in  the  gross 
features,  though  perhaps  not  in  the  details  and  in  the  overall  flux  level.  The  theoretical 
spectrum  calculated  by  Habing  (1968)  is  valid  over  the  wavelength  range  912  -  3000  A. 
The  long-dashed  curve  in  Figure  9  shows  his  spectrum  multiplied  by  a  factor  of  1.6,  which 
accounts  for  the  effects  of  a  different  distance  and  attenuation  between  our  geometry  and 
his.  Draine’s  (1978)  spectrum  for  the  wavelength  band  of  912  -  2000  A  is  based  on  the 
compilation  of  observational  and  theoretical  data.  Our  spectrum  agrees  quite  well  with 
his  within  this  wavelength  range  if  we  multiply  his  flux  by  a  factor  of  1.7  (short  dashed 
curve),  again  accounting  for  the  different  geometries  in  the  two  cases. 

5.  The  Interaction  Network 

If  the  dark  matter  in  Draco  is  baryonic — whether  distributed  continuously  through¬ 
out  the  system,  or  composed  of  many  small  cloud  constituents — it  is  subject  to  several 


15 


key  chemical  interactions,  which  we  here  describe.  In  the  interest  of  keeping  our  calcula¬ 
tions  simple  yet  accurate,  we  have  chosen  to  use  elements,  species  and  reactions  based  on 
those  seen  to  dominate  in  typical  interstellar  clouds.  Our  network  is  therefore  primarily 
composed  of  the  reactions  and  elements  used  by  Roberge  (1981).  Even  though  this  net¬ 
work  IS  relatively  simple,  it  produces  results  that  generally  agree  with  the  more  elaborate 
chemical  network  of  Black  and  Dalgarno  (1977).  We  also  use  some  components  of  the 
detailed  chemical  networks  employed  by  VanDishoeck  and  Black  (1986),  and  Tielens  and 
Hollenbach  (1985). 

Our  simple  chemical  network  is  composed  of  the  elements:  hydrogen,  helium,  carbon, 
nitrogen,  oxygen,  iron,  silicon,  sulfur  and  magnesium.  These  elements  are  used  in  70 
different  reactions  to  produce  and  destroy  a  total  of  32  different  species,  which  include: 
H,  He,  C,  N,  O,  Si,  S,  Mg,  Fe,  free  electrons,  H2,  OH,  CO,  O2,  H2O,  H+,  He+,  C+, 
N+,  0+,  Si+,  S+,  Mg+,  Fe+,  H+,  OH+,  CO+,  H2O+,  HCO+,  H+,  H3O+  and  H2*  (i.e., 
the  vibrational-rotationally  excited  ground  electron  state  of  H2).  A  complete  list  of  the 
reactions  employed  in  these  calculations  is  given  in  Table  3. 

The  reactions  of  our  network  can  be  classified  as  positive  ion  (PI);  neutral  (NR);  elec¬ 
tron  recombination  (El);  radiative  electron  recombination  (AI);  cosmic  ray  (CR);  photoion¬ 
ization  and  dissociation  (PP);  hydrogen-grain  (H2)  and  rotational-vibrationally  excited 
hydrogen  (MH)  reactions  with  a  characteristic  associated  rate  specific  for  each  reaction. 
These  rates  are  generally  dependent  on  parameters  such  as  the  cloud  temperature,  FUV 
flux,  line  and  continuum  attenuation  of  the  FUV  flux  and  the  intensity  of  cosmic  rays 
within  the  dwarf  Spheroidal  galaxy.  To  reflect  the  importance  of  each  of  these  parameters 
Table  3  also  lists  the  rate  coefficients  for  each  reaction.  The  rates  per  particle  number 
density  for  the  PI,  NR,  El  and  AI  reactions,  which  are  primarily  governed  by  the  cloud 
temperature,  can  all  be  calculated  with  the  equation 

k  =  A  (r/300  K)“  exp(-/9/T)  cm^  s’^  ,  (7) 

where  T  is  the  temperature  of  the  cloud  in  Kelvins  and  A,  a  and  ^  are  the  rate  coefficients 
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for  each  specific  reaction.  As  for  the  CR,  reactions,  their  rate  is  given  by 

k  =  CA,  (8) 

where  A  is  the  rate  coefficient  for  each  specific  reaction  and  =  1.3  x  10“^^  s~^  for 
conditions  in  the  interstellar  medium.  We  expect  to  be  much  smaller  within  Draco. 
The  rates  for  the  PP  reactions  are  based  on  a  plane  parallel  assumption  for  the  radiation 
transfer  and  an  interstellar  FUV  flux  following  Habing  (1968): 

k  =  AS  flux  exp{-^Au)  s~^  ,  (9) 

where  Ay  is  the  visual  extinction,  which  increases  with  depth  into  the  cloud,  and  is  pri¬ 
marily  due  to  dust  or  other  continuum  extinction  processes.  A  and  ^  are  rate  coefficients 
for  each  specific  reaction,  and  the  parameter  S fiux  is  a  factor  that  scales  the  level  of  FUV 
flux  within  Draco  relative  to  Habing’s  results  for  the  interstellar  medium  (see  §  4  above). 
This  simple  scaling  is  possible  because  the  rates  are  directly  proportional  to  this  flux. 

The  rate  coefficients.  A,  a  and  ^  in  each  of  the  above  rate  equations  were  derived  by 
Prasad  &  Huntress  (1980)  based  on  theoretical  and  empirical  results  and  are  included  in 
Table  3  for  every  reaction.  This  dependence  on  the  extinction  is  based  on  empirical  fits 
to  observations  along  lines  of  sight  within  the  galactic  plane.  If  the  typical  grain  size  in 
Draco  is  bigger  than  that  within  our  Galaxy,  then  the  extinction  law  must  be  modified.  It 
is  beyond  the  scope  of  the  work  reported  here  to  reproduce  in  detail  an  attenuation  curve 
for  a  variety  of  dust  grain  distributions  analogous  to  that  shown  for  the  ISM  in  Figure  6. 
Instead,  we  shall  try  to  obtain  qualitative  results  by  assuming  that  the  attenuation  may 
still  be  written  in  terms  of  Ay  with  the  rate  coefficients  as  shown,  though  the  dependence 
of  Ay  on  the  dust  constituents  is  modified  by  the  different  grain  radius  r ^ .  When  compared 
with  the  visual  extinction  in  the  Galaxy,  the  value  of  Ay  in  Draco  will  be  enhanced  by  a 
factor  (rg/0.5  ^m)^  corresponding  to  the  grain  cross  section,  and  a  factor  (0.5  ixm/rgY 
for  the  grain  number  density,  under  the  assumption  that  the  total  amount  of  C,  0,  and  N 
content  of  the  grains  is  fixed  to  be  about  half  of  the  observed  abundance  in  this  system 
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and  that  most  of  the  heavier  elements  are  contained  therein  (see  §  2  above).  For  simphcity, 
we  shall  also  assume  that  the  optical  efficiency  in  the  wavelength  range  of  interest  scales 
approximately  as  {rg/0.5  //m)^.  This  treatement  of  the  dust  attenuation  will  be  improved 
in  future  work. 

In  this  network,  we  concentrate  our  efforts  on  the  species  of  H,  H2  and  H2  since  these 
are  the  most  abundant  and  therefore  the  most  likely  to  be  detectable  at  the  distance  of 
Draco.  This  entails  a  detailed  investigation  of  the  formation  and  destruction  of  H,  H2  and 
H2  and  their  effect  on  other  species  within  the  cloud.  One  primary  process  for  forming 
H2  is  the  absorption  of  a  Lyman  or  Werner  band  photon  which  initially  excites  H2  from 
its  ground  electron  state,  X^S+,  to  a  higher  electron  state  of  (for  Lyman)  or 

(for  Werner),  followed  by  a  fluorescence  to  a  bound  vibrational-rotationally  excited  ground 
electron  state,  H2  (Black  &  VanDishoeck  1987;  Martin,  Hurwitz  &  Bowyer  1990).  This 
reaction  is  listed  as  item  MH14  in  Table  3E  and  its  rate  is  given  by  Tielens  &  Hollenbach 
(1985)  as 

k  =  3.06  X  10-'«  Sfiux  l^{r)  exp(-2.5  A,)  s"'  .  (10) 

This  rate  takes  into  account  the  continuum  attenuation  due  to  dust  through  the  visual 
extinction  Ay  which,  again,  increases  with  depth  into  the  cloud  and  is  modified  from  its 
galactic  value  as  described  above.  Using  the  factor  this  rate  equation  also  takes  into 
account  the  H2  line  self-shielding  of  the  Lyman  and  Werner  band  photons.  When  photon 
absorption  is  dominated  by  the  Doppler  cores  or  the  Lorentz  wings  (r  >  10),  the  H2  line 
self-shielding  factor  is  given  by 

/9(r)  =  {r“^  [ln(T/7r^/2)]~^/2  ^  (b/T)^/^}  eT{c{TbTV~'^v^'^Y^‘^  ,  (11) 

where  the  dimensionless  parameters  b,  Vi  and  r  are  given  by  6  =  9.2  x  10“^^t;j\  vi  = 
5  X  10  ,  r  =  1.2  X  10  N [H2)6v^  ^ ,  and  equals  the  turbulent  Doppler  line  width 

in  km  s  and  N(^H2)  is  the  column  density  of  H2.  For  the  linear  part  of  the  curve  of 
growth  (r  <  10),  the  H2  line  self- shielding  factor  is  given  by 

/3(r)  =  S(-l)”  T”/|n!(n  +  W"]  ,  (12) 
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where  the  summation  is  over  n  from  0  to  oo.  Since  ^(r)  is  inversely  proportional  to  r, 
which  is  directly  proportional  to  N(H2),  it  decreases  with  depth  into  the  cloud.  As  with 
the  PP  reactions,  Sfiux  is  the  factor  that  scales  the  FUV  flux  relative  to  that  used  by 
Habing  (1968). 

Black  &  VanDishoeck  (1987)  have  shown  that  another  source  for  producing  is 
the  merging  of  two  hydrogen  atoms  on  a  dust  grain.  When  the  grain  number  density 
IS  1  X  10“^2  times  the  number  density  of  hydrogen  nuclei  (n^)— a  typical  value  in  the 
interstellar  medium — the  rate  for  this  reaction  is 


k  =  3x  10-18T1/2  yf  uh  n(H)  cm-^  s"'  , 


(13) 


where  yf  is  the  efficiency  parameter  (usually  3)  and  n(H)  is  the  number  density  of  H. 

The  destruction  of  H2  is  primarily  due  to  de-excitation  by  the  spontaneous  emission 
of  infrared  radiation,  collisions  and  reactions  with  other  species,  and  the  absorption  of  a 
FUV  photon.  The  rates  per  particle  number  density  for  most  of  these  reactions  are  based 
on  the  rate  coefficients  reported  by  Tielens  &  Hollenbach  (1985)  (see  Table  3),  and  are 
given  by 


k  =  A(T/300  K)“  exp{-/3/T)  cm^  s"'  .  (14) 

The  peculiar  rate  of  reaction  MH3  (involving  and  H2)  and  that  of  reaction  MHl  (the 
absorption  of  a  FUV  photon)  are  also  given  by  Tielens  &  Hollenbach  (1985)  as 


k  =  lAx  10-12  yi/2  exp[-18100/(r  +  1200)]  cm^  s”!  ,  (15) 

and 

k  =  1.0  X  10-11  Sfiux  exp(-2.5  A„)  s"!  ,  (16) 

respectively. 

H2  IS  formed  by  many  process,  as  is  evident  from  the  network  of  reactions  listed  in  Table 
3,  but  it  is  its  destruction  by  a  Lyman  or  Werner  band  photon  that  affects  its  population 
the  most.  And  when  H2  is  a  dominant  species  in  a  cloud,  this  destruction  process  is  the 
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main  source  for  producing  H.  The  rate  for  this  reaction  is  given  by  Tielens  k  Hollenbach 
(1985)  as 

k  =  3Ax  10~^^  S flux  0iT)  exp(-2.5A„)  s“^  ,  (17) 

where  again  the  shielding  factor  ^(r)  is  used. 


6.  Calculational  Procedure  and  Results 


The  dark  matter  distribution  we  derived  in  §  3  above  specifies  the  average  density 
(or  equivalently,  the  enclosed  mass)  as  a  function  of  radius.  Its  dumpiness,  however,  is 
not  known.  We  shall  therefore  consider  two  extreme  cases,  one  in  which  the  matter  is 
distributed  continuously  throughout  the  dSph  (i.e.,  no  dumpiness),  and  a  second  in  which 
the  gas  is  condensed  into  Me  =  10  Mq  clouds.  A  more  general  profile,  constituting  a  range 
of  cloud  sizes  and  masses  will  be  left  to  future  work. 

In  both  cases,  the  cloud  temperature  T  is  assumed  to  be  the  representative  value  of 
100  K  attributed  to  most  known  clouds.  In  the  first  model — that  of  the  continuously 
distributed  gas  the  dark  matter  density  is  simply  given  by  the  best  fit  PK  profile 

used  by  Haller  &  Melia  (1996),  which  is  shown  graphically  in  Figure  2.  A  lower  limit  to  the 
size  of  the  clouds  in  the  second  model  is  the  Jeans  radius  Rj  (Jeans  1902),  below  which 
the  gas  condensation  will  have  collapsed  and  possibly  formed  stars.  In  terms  of  the  cloud 
mass  Me  and  mean  molecular  weight  p,  this  limiting  radius  is  given  as 


„  _  3/i  ruH  GMe 

n^kT  ’ 
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where  ttiu  is  the  proton  mass.  For  simplicity,  we  shall  adopt  this  value  as  the  characteristic 
radius  of  all  the  (spherical)  clouds  in  the  system,  and  assume  that  they  have  uniform 
average  density,  given  by 


<  Pe 


ZMr 


(19) 


Defining  ne{r)  to  be  the  cloud  number  density  as  a  function  of  radius,  it  is  clear  that  in 
this  case 

Pdm  =  ne{r)Me  .  (20) 
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Numerically,  for  10  M©  clouds,  i?j  ~  2  x  10  ^  pc  and  <  >~  3  x  10®  M©  pc~®.  For 

these  cloud  characteristics,  the  volume  filling  factor  in  Draco  would  be  about  2.8  x  10~®, 
with  an  areal  filling  factor  of  roughly  1  x  10"®.  From  an  observational  point  of  view,  only 
about  0.1  %  of  Draco  would  therefore  be  obscured  by  all  the  clouds  viewed  as  a  whole, 
which  would  be  consistent  with  the  observed  lack  of  extinction  within  this  system.  The 
cloud  number  density  is  about  9  x  lO”^  pc“*  at  the  center  of  the  dSph  and  falls  off  rapidly 
towards  zero  as  the  tidal  radius  is  approached.  Thus,  a  typical  spacing  of  the  clouds  in 
the  center  of  Draco  is  approximately  2  pc. 

As  discussed  in  §  2  above,  the  observed  average  dust  grain  size  in  many  Galactic  loca¬ 
tions  appears  to  be  about  0.5  ^m,  but  because  the  important  dust  destruction  mechanisms 
that  limit  the  grain  size  are  presumably  not  functioning  in  dSphs,  we  consider  values  as 
high  as  48  pm  in  our  calculations.  We  also  assume  that  the  cosmic  ray  intensity  in 
Draco  is  much  lower  than  in  the  Galactic  plane,  and  we  take  the  FUV  flux  to  be  about 
1.6  (=  Sfiux)  times  as  large  as  Habing’s  (1968)  value  on  the  basis  of  our  discussion  in  §  4 
above.  Finally,  the  gas  turbulent  velocity  dispersion  is  taken  to  be  the  usual  ~  1  km  s~^ 
typically  seen  in  most  interstellar  clouds  in  the  solar  neighborhood  for  the  clumpy  case, 
and  equal  to  the  total  velocity  dispersion  in  Figure  3  for  the  continuous  case. 

The  calculational  procedure  calls  for  the  simultaneous  solution  of  32  non-linear  differ¬ 
ential  equations  to  determine  the  equilibrium  density  of  each  species  at  discrete  steps  into 
the  cloud(s).  Using  the  interaction  network  developed  in  §  5  above,  all  of  these  equations 
may  be  written  in  the  generic  form  (for  each  species  i) 

^  =  E  +  E  + E  ‘if  "i  -  "i  E 

i  j  j 

x:  if  -  -  n,  Y.  ’’f-  .  (21) 

i  3 

where  superscripts  “-f”  and  stand  for  creation  and  annihilation  processes,  respec¬ 
tively,  gp  IS  for  gas  processes,  pp  is  for  photo-processes,  and  cr  stands  for  cosmic  ray. 
In  equilibrium,  dtijldt  0,  and  the  equations  may  be  solved  simultaneously  using  the 


21 


Newton-Raphson  method.  Starting  with  the  full  (unattenuated)  FUV  flux  at  the  cloud 
surface,  we  integrate  inwards,  taking  into  account  the  gradual  extinction  of  the  radiation, 
until  the  densities  of  the  various  species  reach  their  asymptotic  (i.e.,  as  r  ^  0)  values. 
This  steady  state  is  brought  on  by  the  attenuation  of  the  incoming  radiation  by  the  dust 
and  H2  to  the  point  where  the  reaction  rates  at  smaller  radii  are  affected  primarily  by  the 
internal  energy  of  the  clouds.  In  every  case,  the  radiation  is  attenuated  significantly  and 
the  species  reach  their  asymptotic  values  within  a  small  fractional  percentage  of  the  cloud’s 
radius  near  the  surface,  justifying  a  plane  parallel  treatment  of  the  radiative  transfer  in 
this  region  (see  below). 

In  Figure  10,  we  show  the  equilibrium  distribution  of  the  dominant  species  (H,  H2, 
and  H2)  as  a  function  of  the  radius  (in  pc)  within  Draco,  for  the  case  of  a  continuously 
distributed  gas.  The  FUV  flux  induces  rapid  changes  near  the  dSph’s  surface  until  the 
rapidly  increasing  visual  extinction  into  the  cloud  and  the  H2  self-shielding  quenches  the 
radiation  field  and  the  gas  converts  almost  entirely  to  H2 .  The  assumed  grain  size  for  this 
simulation  is  48  //m.  This  distribution  has  two  major  difficulties:  (1)  the  corresponding 
total  atomic  hydrogen  content  is  Mh  ^  10®  M©,  and  (2)  the  visual  extinction  through 
the  galaxy  would  be  Av  1.7  x  10®.  Both  of  these  are  inconsistent  with  the  observations, 
which  indicate  that  Draco  is  effectively  transparent  and  that  it  contains  at  most  ~  68  M© 
of  atomic  hydrogen  (see  §  2). 

By  comparison,  these  difficulties  are  not  present  when  the  gas  is  clumpy,  say  in  the 
form  of  10  M©  clouds.  The  equilibrium  profiles  for  this  case  are  shown  in  Figures  11a), 
b),  and  c)  assuming  the  same  FUV  flux  and  dust  content  as  those  above.  For  a  hydrogen 
column  density  Njj  to  the  center  of  each  cloud  (though  concentrated  primarily  near  its 
surface),  the  total  hydrogen  content  of  Draco  may  be  estimated  according  to  the  equation 

47r  Rl  Nh  ruH  ,  (22) 

where  Mdm  is  the  total  mass  of  dark  matter  in  this  system,  and  Me  and  Rc  are  the 
individual  cloud  mass  and  radius,  respectively.  For  this  simulation,  we  find  that  Nfj 
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9.7  X  10^  cm  for  which  Mu  rs  67  Mq.  We  note  that  the  corresponding  visual  extinction 
through  each  cloud  is  w  1.1  x  10^.  However,  as  noted  earlier,  the  areal  filling  factor  for 
this  case  is  only  ~  0.1%  so  that  typical  lines-of-sight  through  this  system  would  presumably 
not  encounter  any  significant  extinction. 

Reducing  the  dust  grain  size  drives  the  equilibrium  away  from  H2  towards  H  since  the 
attenuation  of  the  external  FUV  radiation  is  then  not  as  effective.  So,  for  example,  if  we 
reduce  rg  to  ~  0.7  ^m,  Draco  would  then  need  about  500  times  more  dust  in  order  not 
violate  the  21  cm  limit  on  H.  We  will  discuss  the  various  observational  implications  of 
these  results  in  the  next  section. 

7.  Observational  Implications 

On  the  basis  of  the  21  cm  data,  it  would  appear  that  the  dark  matter  in  Draco  could 
be  mostly  H2,  but  only  if  it  is  in  the  form  of  ~  10  —  100  Mq  clouds  and  if  the  dust 
grains  are  typically  at  least  ~  90  times  larger  than  those  in  the  Galaxy.  This  might  be 
consistent  with  the  fact  that  the  dominant  size-dependent  dust  destroying  mechanisms  in 
our  galaxy  are  absent  in  typical  dSphs.  Interestingly,  the  absence  of  dynamical  processes 
(such  as  supernova  explosions)  that  would  normally  induce  massive  star  formation  (which 
is  non-existent  in  Draco)  may  also  account  for  this  difference  in  the  dust  structure.  A 
continuously  distributed  gas  is  ruled  out  by  the  lack  of  any  observed  extinction  in  this 
system,  which  would,  however,  be  consistent  with  the  very  low  areal  covering  fraction 
(~  0.1  -  1%)  if  the  dust  and  gas  were  concentrated  in  clouds. 

We  may  estimate  a  gross  upper  limit  to  the  FIR  flux  produced  by  the  dust  by  simply 
taking  the  intercepted  UV  radiation  and  assuming  that  it  is  completely  reprocessed  into 
this  longer  wavelength  band.  Thus,  using  a  FUV  flux  of  Ri  5.3  x  10“^  ergs  cm“2  at 
Draco,  the  k,  10®  10  M©  clouds,  each  with  a  radius  of  «  2  x  10~^  pc,  would  intercept 
a  total  radiative  power  of  6  x  10^'^  ergs  s“^  At  earth,  a  conservative  upper  limit  to  the 
FIR  flux  from  dust  would  therefore  be  1  x  10”^®  ergs  cm~^  or  approximately  0.3  Jy 
if  it  is  all  emitted  near  10  ^m.  By  comparison,  the  IRAS  sensitivity  at  this  wavelength 
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is  w  0.07  Jy  (Fazio  k  Eisenhardt  1990).  In  reality,  only  a  small  fraction  of  the  FUV  flux 
can  actually  be  reprocessed  in  this  fashion,  and  so  this  number  must  be  a  very  generous 
estimate.  That  is,  observations  with  current  instrimients  probably  would  not  be  able  to 
detect  this  dust  emission.  SIRTF  is  projected  to  have  a  sensitivity  of  6  x  10“®  Jy  at  10  yum, 
but  it  too  may  not  be  able  to  detect  this  dust  since  its  field  of  view  is  too  small  to  observe 
more  than  a  small  fraction  of  the  entire  dSph. 

The  situation  also  does  not  appear  to  be  promising  when  it  comes  to  the  H2  emission 
within  the  J  and  K  bands.  In  Figures  12a)  and  b)  we  show  the  expected  spectrum  (K 
and  J  infrared  bands,  respectively)  integrated  over  the  entire  volume  in  Draco,  as  seen  at 
earth.  These  are  based  on  the  236  strongest  infrared  emission  lines  of  excited  H2  ranging 
in  wavelength  from  12  fim  to  0.85  ^m,  as  described  by  Black  &  VanDishoeck  (1987). 
NICMOS  on  HST  will  be  an  assembly  of  3  cameras,  with  various  fields  of  view  (FOV) 
and  a  range  of  filters  (Thompson  1996).  Our  predicted  line  flux  in  the  K  band  is  typically 
~  2  X  10"^^  ergs  cm“2  g-i,  but  with  its  actual  FOV,  Camera  2  (which  is  sensitive  in  the 
wavelength  range  1.75  -  2.35  /urn)  will  be  able  to  sample  about  10~®  of  Draco’s  volume 
with  a  pointed  integration  at  its  center.  At  this  level,  it  would  not  be  able  to  detect  the 
clouds  in  any  reasonable  integration  time. 

In  the  J  band,  the  expected  line  flux  for  the  most  prominent  components  is  typically 
~  9  X  10  ergs  cm~^  s~^,  and  using  the  same  assumptions  for  Camera  1  (which  is 
sensitive  in  the  wavelength  range  1.0  -  1.8  //m,  has  a  slightly  smaller  FOV  and  should 
be  able  to  sample  roughly  8  x  10“®  of  Draco’s  volume),  none  of  these  lines  will  be  easily 
detectable.  The  prospects  for  detection  of  these  J  band  lines  are  similar  with  Camera  3, 
even  though  it  has  a  larger  FOV. 

8.  Concluding  Remarks 

Interestingly,  similar  considerations  as  the  ones  we  have  described  here  were  sketched 
out  earlier  by  Jura  (1977)  for  elliptical  galaxies.  He  speculated  that  because  the  UV 
radiation  in  these  systems  is  low,  the  interstellar  clouds  in  ellipticals  may  be  cold  and 
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largely  H2.  As  a  result,  the  very  low  mass  clouds  (A/  ^  10  Mq)  may  collapse  to  create 
stars,  and  only  late  type  stars  may  form.  The  conditions  we  have  explored  for  the  dSph 
galaxies  are  not  unlike  these.  Our  conclusions  would  seem  to  indicate  that  star  formation 
within  systems  such  as  Draco  may  be  largely  dependent  on  which  clouds  can  collapse 
without  substantial  external  influence.  If  supernova-induced  shocks  or  large-scale  density 
waves  are  instrumental  in  producing  large  mass  condensations  that  lead  to  the  formation 
of  early  type  stars,  the  general  absence  of  these  in  dSphs  may  be  consistent  with  the  fact 
most  of  the  clouds  in  these  galaxies  may  not  have  undergone  this  type  of  dynamically- 
induced  collapse.  As  a  result,  most  of  the  mass  (perhaps  as  much  as  99%)  may  still  be  in 
the  form  of  primordial  H2  clouds  and  only  late  type  stars  formed  from  the  few  lower  mass 
clouds  that  crossed  below  the  Jeans  criterion. 

Finally,  we  check  to  see  whether  such  a  “cluster”  of  H2  clouds  could  in  fact  have  avoided 
core  collapse  over  a  Hubble  time  tn-  Since  the  clouds  presumably  dominate  the  mass  of 
the  system,  we  may  use  the  simple  estimate  (Binney  k  Tremaine  1987)  of  the  core  collapse 

time  scale  tec  ~  12  —  19  for  a  single  mass  component  in  terms  of  the  median  relaxation 
time 

,  O.UNc 

trh  =  - 

ln(0.4iVe) 

where  is  the  number  of  clouds,  r/,  is  the  system’s  median  radius  and  Mdm  is  the  total 
mass  of  dark  matter.  For  Draco,  tec  ~  1  X  10^^  —  1.5  x  10^^  years,  much  longer  than  tn- 
This  research  was  supported  by  NSF  (PHY  88-57218)  and  NASA  (NAGW-2822).  We 
are  very  grateful  to  R.  Kuruez  for  sending  us  his  calculated  stellar  spectra,  and  to  E. 
Olszewski  for  informative  discussions  and  for  providing  the  velocity  data  for  Draco. 
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Luminosity  Class  V  -  YII  (main  sequence,  sub  dwarfs  and  white  dwarfs) 

model  #  Spectral  Type  TefF  log(g)  log(Rs/R@)  [M/H]*  Kurucz  model  #  Kurucz  TefF  Kurucz  lo] 
D1  05(05-09)  47000  4  1.25  0  408  47500  5 
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Luminosity  Class  I  -  II  (super  giants  and  bright  giants) 

model#  Spectral  Type  Teff  IpgCg)  Iog(Rs/Ho)  [M/H]*  Kurucz  model#  Kurucz'^efF  Kuruczloi 
si  05(05-09)  -42000  3.4*  -1.35  0  406  ^  42500  5~ 

S2  B0(B0-B4)  25500  3.1  1.3  Q  365  26000  3 
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TABLE  2A 


Luminosity  Class  V  -  VII  (main  seguence,  sub  dwarfe  and  white  dwarfs) 
model  #  Spectral  Type _ percentage _ number  of  stars* 


D1 

D5(0§-G9) 

-0.4)00048^35 

4.89350E+04 

D2 

B0(B0-B4) 

0.012532 

1.25320E+07 

D3 

B5(B5-B9) 

5.0004 

5.00040E+09 

D4 

A0(A0-A4) 

5.4024 

5.40240E+09 

D5 

A5(A5-A9) 

2.3153 

2.31530E+09 

D6 

F0(F0-F4) 

3.4817 

3.48170E+09 

D7 

F5(F5-F9) 

3.4817 

3.48170E+09 

D8 

G0(G0-G4) 

4.444 

4.44400E+09 

D9 

G5(G5-G9) 

4.444 

4.44400E+09 

DIO 

K0(K0-K4) 

4.8373 

4.83730E+09 

Dll 

K5(K5-K9) 

4.8373 

4.83730E+09 

D12 

M0(M0-M4) 

30.376 

3.03760E+10 

D13 

M5(M5-M9) 

30.376 

3.03760E+10 

*Based  on  a  total  of  lE+1 1  stars  in  the  galaxy. 
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TABLE  2B 


Luminosity  Class  III  -  IV  (giants  and  subgiants) 


model  #  Spectral  Type _ percentage _ number  of  stars* 


G1 

05(05-09) 

0.000000434 

4.34000E+02 

G2 

B0(B0-B4) 

2.591 7E-06 

2.591 70E+03 

G3 

B5(B5-B9) 

0.0010341 

1.03410E+06 

G4 

A0(A0-A4) 

0.016879 

1.68790E+07 

G5 

A5(A5-A9) 

0.0072339 

7.23390E+06 

G6 

F0(F0-F4) 

0.10455 

1.04550E+08 

G7 

F5(F5-F9) 

0.10455 

1.04550E+08 

G8 

G0(G0-G4) 

0.12361 

1.23610E+08 

G9 

G5(G5-G9) 

0.12261 

1.22610E-K)8 

GIO 

K0(K0-K4) 

0.23635 

2.36350E+08 

Gil 

K5(K5-K9) 

0.23635 

2.36350E+08 

G12 

M0(M0-M4) 

0.01742 

1.74200E+07 

G13 

M5(M5-M9) 

0.01742 

1.74200E+07 

*Based  on  a  total  of  lE+1 1  stars  in  the  galaxy. 


32 


TABLE  2C 


Luminosity  Class  I  -  II  (super  giants  and  bright  giants) 


model  #  Spectral  Type _ percentage _ number  of  stars* 


SI 

05(05-09) 

0.00000027 

2.70000E+02 

S2 

B0(B0-B4) 

1.2565E-06 

1.25650E+03 

S3 

B5(B5-B9) 

3.7695E-06 

3.76950E+03 

S4 

A0(A0-A4) 

0.000073819 

7.38190E+04 

S5 

A5(A5-A9) 

0.000031637 

3.16370E+04 

S6 

F0(F0-F4) 

0.00022851 

2.285  lOE+05 

S7 

F5(F5-F9) 

0.00022851 

2.285  lOE+05 

S8 

G0(G0-G4) 

0.00032297 

3.22970E+05 

S9 

G5(G5-G9) 

0.00032297 

3.22970E+05 

SIO 

K0(K0-K4) 

0.00056577 

5.65770E+05 

Sll 

K5(K5-K9) 

0.00056577 

5.65770E+05 

S12 

M0(M0-M4) 

0.0003636 

3.63600E+05 

S13 

M5(M5-M9) 

0.0003636 

3.63600E+05 

*Based  on  a  total  of  lE+1 1  stars  in  the  galaxy. 
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TABLE  3A 


GAS  PHASE  REACTIONS 


Number 

Reaction 

A 

(X 

3 

PI15 

He+  +  H2->H+  +  H  +  He 

1.5(-13) 

0.00 

- K - 

0.00 

PI49 

H2+  +  H  ^  H+  +  H2 

I.O(-IO) 

0.00 

0.00 

PI43 

0+  +  H2  ^  0H+  +  H 

1.6(-09) 

0.00 

0.00 

PI42 

0+  +  H  ^  0  +  H+ 

7.0(-10) 

0.00 

0.00 

PI13 

He+  +  H-4-He  +  H+  +  hv 

1.9(-15) 

0.00 

0.00 

PI12 

H+  +  H20^H  +  H20+ 

8.2(-09) 

0.00 

0.00 

PIl 

H+  +  0^H  +  0+ 

7.0(-10) 

0.00 

232.0 

PI32* 

C+  +  OH  ->  CO  +  H+ 

7.7(-10) 

0.00 

0.00 

PI32 

C+  +  OH  ^  C0+  +  H 

7.7(-10) 

0.00 

0.00 

PI50 

H2+  +  H2-4>H3+  +  H 

2.1  (-09) 

0.00 

0.00 

PI19 

He+  +  C0->He  +  C+  +  0 

1.6(-09) 

0.00 

0.00 

PI65 

H3+  0  ->  0H+  +  H2 

8.0(-10) 

0.00 

0.00 

PI59 

0H+  +  H2  ->  H2O+  +  H 

1.1  (-09) 

0.00 

0.00 

PI79 

H2O+  +  H2  ^  H3O+  +  H 

6.1(-10) 

0.00 

0.00 

PI38 

C+  +  H2O  ->  HCO+  +  H 

2.7(-09) 

0.00 

0.00 

PI61 

C0+  +  Hj  ^  HCO+  +  H 

2.0(-09) 

0.00 

0.00 

PI69 

H3+  +  CO  ^  HCO+  +  H2 

1.7(-09) 

0.00 

0.00 

PI27 

C+  +  O2  -4-  CO+  +  0 

7.5(-10) 

0.00 

0.00 

PI28 

C+  +  O2  ^  CO  +  0+ 

4.1(-10) 

0.00 

0.00 

PI16 

He+  +  O2  He  +  0  +  0+ 

1.0(-09) 

0.00 

0.00 

NR4 

0  +  OH  ^  O2  +  H 

3.0(-ll) 

-.36 

0.00 

Rate  coefficients  are  from  Prasad  &  Huntress  (1980),  unless  otherwise  stated. 
Number  notation:  1.0(-10)  =  l.OxlO-i® 

*From  Roberge  (1981). 
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TABLE  3B 


ELECTRON  RECOMBINATION  REArTTONS 


Number 

_ Reaction _ 

A 

(X 

R 

All 

H+  +  e'  H  +  hv 

3.5(-12) 

-.70 

- M - 

0.00 

EI12 

H3+  +  e'  ->  Hj  +  H 

1.1  (-07) 

-.50 

0.00 

EI13 

H3+  +  e‘  H  +  H  +  H 

1.1  (-07) 

-.50 

0.00 

AI5 

O'*'  +  e"  0  +  hv 

3.4(-12) 

-.63 

0.00 

EI7 

OH+  +  e‘  ^  0  +  H 

2.0(-07) 

-.50 

0.00 

EI18 

H2O+  +  e'  ^  OH  +  H 

2.0(-07) 

-.50 

0.00 

EI38 

H3O+  +  e"  ^  OH  +  H  +H 

6.5(-07) 

-.50 

0.00 

EI37 

H3O+  +  e"  ^  H2O  +  H 

6.5(-07) 

-.50 

0.00 

AI3 

C^  +  e  ^  C  +  hv 

4.4(-12) 

-.61 

0.00 

AI2 

He+  +  e-  ^  He  +  hv 

4.5(-12) 

-.67 

0.00 

AI4 

N"*"  +  e"  N  +  hv 

3.8(-12) 

-.62 

0.00 

AI9 

Mg'*'  +  e"  Mg  +  hv 

2.8(-12) 

-.86 

0.00 

AI6 

Si'*'  +  e"  — >  Si  +  hv 

4.9(-12) 

-.60 

0.00 

AI7 

S'*"  +  e"  ->  S  +  hv 

3.9(-12) 

-.63 

0.00 

AI8 

Fe'*’  +  e"  ->  Fe  +  hv 

3.7(-ll) 

-.65 

0.00 

EI8 

CO+  +  e‘  C  +  0 

1.8(-07) 

-.50 

0.00 

EI20 

HCO+  +  e‘^H  +  CO 

2.0(-07) 

-.75 

0.00 

Rate  coefficients  are  from  Prasad  &  Huntress  (1980),  unless  otherwise  stated. 
Number  notation:  L0(-10)=  l.OxlO-io 
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TABLE  3C 


COSMIC  RAY  REACTIONS 


Number 

_ Reaction _ 

A 

a 

B 

CRl 

H  +  cr  H+  +  e- 

4.6(-01) 

0.00 

- - 

0.00 

CR2 

He  +  cr  ^  He'*'  +  e" 

5.0(-01) 

0.00 

0.00 

CR7 

H2  cr  — ^  H2^  0 

9.3(-01) 

0.00 

0.00 

CR6 

Hj  +  cr  ->  H  +  H"'"  +  e' 

2.2(-02) 

0.00 

0.00 

CR8 

Hj  +  cr^H  +  H 

l.O(-Ol) 

0.00 

0.00 

CR3 

C  +  cr  ^  C^  +  e“ 

1.8(-00) 

0.00 

0.00 

CR5 

0  +  cr  ^  0+  +  e‘ 

2.8(-00) 

0.00 

0.00 

CRIO 

CO  +  cr  C0+  +  e' 

3.0(-00) 

0.00 

0.00 

CR4 

N  +  cr  ->  N'*'  +  e" 

2.1  (-00) 

0.00 

0.00 

Rate  coefficients  are  from  Prasad  &  Huntress  (1980),  unless  otherwise  stated. 
Number  notation:  1.0(-10)  =  l.OxlO-io 
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TABLE  3D 


mOTOIONIZATION  AND  DISSOCTATION  REACTTONS 


Number 

_ Reaction _ 

A 

a 

B 

PP55 

+  hv  H  +  H 

(See  text) 

- K - 

PP15 

OH  +  hv  ->  0  +  H 

2.2(-10) 

0.00 

2.0 

PP27 

HjO  +  hv  ^  H2O+  +  e’ 

2.1(-11) 

0.00 

3.1 

PP26 

HjO  +  hv^OH  +  H 

5.1(-10) 

0.00 

1.8 

PPl 

C  +  hv  ->  C*"  +  e' 

2.1(-10) 

0.00 

2.6 

PP17 

CO  +  hv  ^  C  +  0 

1.4(-11) 

0.00 

3.2 

PP5 

Mg  +  hv  ->  Mg'*"  +  e' 

4.4(-ll) 

0.00 

1.4 

PP2 

Si  +  hv  ^  Si"''  +  e‘ 

1.2(-09) 

0.00 

1.6 

PP3 

S  +  hv  +  e~ 

3.9(-10) 

0.00 

2.5 

PP4 

Fe  +  hv  ->•  Fe"*"  +  e' 

I.O(-IO) 

0.00 

2.3 

PP9 

Oj  +  hv  ->  0  +  0 

3.3(-10) 

0.00 

1.4 

Rate  coefficients  are  fi-om  Prasad  &  Huntress  (1980),  unless  otherwise  stated. 
Number  notation:  I.O(-IO)  =  1,0x1 0-io 
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TABLE  3E 


ROTATION-VIBRATIONALLY  EXCITED  HYDROGEN  REACTIONS 


Number 

_ Reaction _ 

A 

a 

3 

H21 

H  +  H  +  Grain  Hj*  +  Grain 

(See  text) 

- - 

MH14 

+  hv  ^  Hj* 

(See  text) 

MHl 

H^  +  hv^H  +  H 

(See  text) 

MH13 

Hj*  ^  +  hv 

2.0(-07) 

0.00 

0.00 

MH2 

Hj*  +  H  ^  Hj  +  H 

1.9(-12) 

.50 

1.0(03) 

MH3 

(See  text) 

MH4 

H2*+H->H  +  H  +  H 

9.8(-12) 

.50 

2.7(04) 

MH5 

9.8(-12) 

.50 

2.7(04) 

MH6 

H2*+H2*  ^H  +  H  +  H2 

9.8(-12) 

.50 

0.00 

MH7 

H2*  +  0  ^  OH  +  H 

9.0(-12) 

1.00 

0.00 

MH8 

H2*  +  OH  ^  H2O  +  H 

3.6(-ll) 

0.00 

0.00 

MH9 

H2*  +  O2  ^  OH  +  OH 

I.O(-IO) 

0.00 

0.00 

Rate  coefficients  are  from  Tielens  &  Hollenbach  (1980),  unless  otherwise  stated. 
Number  notation:  I.O(-IO)  =  l.OxlO-i® 
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FIGURE  CAPTIONS 


Fig.  1.  -  The  aperture  PME  applied  to  stellar  radial  velocity  data  (•)  from  the  dwarf 
spheroidal  galaxy  Draco  (Olszewski  1994).  The  data  points  with  an  open  circle  (o)  include 
an  additional  high  velocity  star  in  the  aperture  PME.  The  data  points  are  not  independent; 
at  a  given  projected  radius,  each  point  contains  all  of  the  stars  within  that  radius.  The 
solid  line  shows  the  computed  aperture  PME  curve  corresponding  to  a  King  model  fit  to 
the  central  velocity  dispersion  and  stellar  surface  density  from  Pryor  &  Kormendy  (1990; 
PK)  with  M/L  =  100.  The  dashed  and  dotted  lines  correspond  to  the  allowed  range  in  the 
central  black  hole  mass,  (1.5-2.0)xl0^  for  a  model  proposed  by  Strobel  &  Lake  (1994; 
SL).  On  the  whole,  the  data  are  more  consistent  with  the  King  model.  (From  Haller  & 
Melia  1996) 

Fig.  2.  -  Dark  mass  density  distribution  for  the  best  fit  PK  model  (see  text  and  Fig. 

!)• 

Fig.  3.  -  Velocity  as  a  function  of  radius  within  Draco:  total  velocity  dispersion 
(solid),  tangential  velocity  dispersion  (short  dash),  radial  velocity  dispersion  (dotted), 
Keplerian  circular  velocity  (long  dash),  and  local  escape  velocity  (dot-dash).  The  vertical 
dashed  line  corresponds  to  the  core  radius  at  117  pc. 

Fig.  4.  -  Schematic  diagram  illustrating  the  simple  model  we  adopt  for  the  Milky 
Way  galaxy.  The  stars  are  assumed  to  lie  uniformly  in  a  thin  disk  (the  ‘stellar’  disk)  with 
a  radius  of  15  kpc,  embedded  within  a  uniform  disk  of  dust  with  a  thickness  of  300  pc. 
Draco  lies  at  a  heliocentric  distance  of  67  kpc,  with  a  galactic  latitude  of  35  degrees  and  a 
galactic  longitude  of  86  degrees. 

Fig.  5.  -  The  percentage  of  galactic  stars  within  each  spectral  type:  main  sequence 
(solid),  giants  (dotted),  supergiants  (dashed). 

Fig.  6.  -  The  average  interstellar  attenuation  curve  of  Savage  &  Mathis  (1979). 
Here,  R\  is  the  ratio  of  the  magnitude  of  attenuation  A.\  to  E(^B  —  V),  as  a  function  of 
wavelength. 

^*8*  ~  The  total  FUV  spectrum  at  Draco  due  to  the  contributions  from  all  the 

stars  in  the  Milky  Way  galaxy,  taking  into  account  the  proper  attenuation  through  the 
galactic  plane. 

Fig.  8.  -  Same  as  Figure  7,  except  now  in  units  of  photons  cm“^  s~^  A  \ 

Fig.  9.  Same  as  Figure  7,  except  now  also  showing  for  comparison  the  theoretical 
spectrum  calculated  for  the  interstellar  medium  by  Habing  (1968)  in  the  wavelength  rcinge 
912  -  3000  A  (long  dashed  curve),  here  multiplied  by  a  factor  of  1.6,  and  Draine’s  (1978) 
spectrum  in  the  wavelength  band  912  —  2000  A,  which  is  based  on  the  compilation  of  obser¬ 
vational  and  theoretical  data  (short  dashed  curve).  Draine’s  spectrum  is  here  multiplied 
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by  a  factor  of  1.7.  In  both  cases,  the  multiplicative  factor  accounts  for  the  dilferences  in 
geometry  (distance  and  degree  of  attenuation)  between  their  case  and  ours. 

Fig.  10.  -  Equilibrium  distribution  of  the  dominant  species  (as  labeled)  as  a  function 
of  the  radius  (in  pc)  within  Draco,  for  the  case  of  a  continuously  distributed  gas.  The 
rapid  changes  near  the  dSph’s  outer  surface  are  due  to  the  effects  of  the  external  FUV 
radiation  field  (mostly  from  the  Milky  Way  galaxy).  The  vertical  dashed  line  indicates  the 
core  radius  at  117  pc. 

Fig.  11.  -  Same  as  Fig.  10,  but  now  for  a  10  M©  cloud  as  described  in  §  6  of  the 
text.  The  three  different  parts  of  this  figure  pertain  to  the  hydrogen-based  species  (a),  the 
oxyS^ii'b^sed  species  (b),  and  the  carbon-based  species  (c). 

Fig.  12.  -  The  predicted  H2  spectrum  integrated  over  the  entire  volume  of  Draco  for 
a)  K  band  and  b)  J  band. 
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Thesis  Abstract,  References  and  Data 


ABSTRACT 


High-precision  velocity  measurements  of  faint  individual  stars  in  the  Galactic  dwarf 
spheroidal  galaxy  Draco  indicate  that  the  mass-to-light  ratio  in  this  system  may  be  as 
high  as  ~  100.  Various  possible  forms  of  dark  matter  have  been  proposed  for  this  class  of 
objects,  but  a  more  recent  analysis  based  on  the  Projected  Mass  Estimator  technique  seems 
to  suggest  that  the  non-luminous  mass  is  distributed  and  may  trace  the  stellar  population. 
Here,  we  examine  in  detail  the  possible  contribution  of  H2  to  the  dark  matter  content  of 
Draco  and  conclude  that  most  of  its  mass  may  be  in  the  form  of  ~  10  —  100  M©  H2  clouds 
if  the  typical  dust  grain  size  is  at  least  ~  90  times  larger  (i.e.,  ^  48  //m)  than  in  our  galaxy. 
This  may  be  consistent  with  the  fact  that  the  dominant  dust  destroying  mechanisms  (e.g., 
supernova  shocks)  that  suppress  the  average  grain  size  in  the  galactic  plane  are  absent  in 
typical  dwarf  Spheroidals.  Thus  the  absence  of  dynamical  processes  (such  as  supernova 
explosions)  that  would  normally  induce  massive  star  formation  (which  is  non-existent  in 
Draco)  may  also  account  for  the  large  grain  size.  The  lack  of  any  observed  extinction  in 
this  system  would  be  consistent  with  the  very  low  areal  covering  fraction  (~  0.1  —  1%) 
of  the  clouds.  We  show  that  the  predicted  H  content  in  Draco  is  then  below  current  21 
cm  detection  limits,  and  that  even  NICMOS  on  HST  would  be  hard  pressed  to  detect  the 
prominent  J  and  K  band  emission  lines  of  excited  H2. 
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